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Abstract 



Abstract 

Gallium Nitride (GaN) is used as an absorbent layer in multijunction cells. These 
solar cells are mainly used in space because the yhave agreat resistance to cosmic 
radiation high-energy (MeV) thanks to the direct and wide energy gap. The only 
drawbackinthis solar cell isweakphoto-current, which affects the overall current of the 
multijunction solar cell. In this work we study the effect of increasing the proportion 
ofindium(In) which replaces Ga on the solar cell properties. The goal is to improve the 
photocurrent to be used in multijunction solar cells based on p-In x Gai. x N/n-In x Gai. x N, 
junction. It was found thatthe proportion ofindium(x = 0.78)greatercurrent(13.51m 
A). However other parameters are degraded. 

Keywords: Gallium Nitride, solar cell (pn), simulation (Silvaco/ Atlas). 
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Resume 

Le nitrure de gallium (GaN) est utilise comme couche absorbante dans des cellules a 
jonctions multiples. Ces cellules solaires sont principalement utilisees dans l'espace, 
car ils ont une grande resistance au rayonnement cosmique de haute energie (MeV) 
grace au gap d'energie direct et large. Le seul inconvenient de cette cellule solaire est 
un faible photo-courant, ce qui affecte le courant global de la cellule solaire a 
jonctions multiples. Dans ce travail, nous etudions l'effet de l'augmentation de la 
proportion d'indium (In) qui remplace le Ga sur les proprietes des cellules solaires. 
L'objectif est d'ameliorer le photocourant pour etre utilise dans des cellules solaires a 
jonctions multiples a base de la jonction p-In x Gai_ x N/n- In x Gai_ x N. On a trouve que la 
proportion d'indium (x = 0,78) donne le maximum du photo-courant (13.51 m A). 
Cependant, les autres parametres sont degrades. 

Mots cles:Nitrure deGallium, cellules solaires (pn), simulation (Atlas/Silvaco). 



ii 



SUMMARY 



Summary : 

Acknowledgments 

Abstract 

Summary 

General Introduction 
Chapter I: Properties of Gallium Nitride 

1.1 Introduction 

1.2 Crystallographic Properties 

1.2.2 The structure of zincblende type cubic: 

1.2.3 The structure NaCl type 

1.3 Chemical properties of GaN 

1.3.1 The surface chemistry of GaN 

1.3.2 The stability of GaN at high temperature 

1.4 Electrical properties of GaN 

1.4.1 GaN mobilities 

1.5 Optical properties of GaN 

1.5 .1 Band gap energy of GaN 

1.5.2 The temperature and pressure dependence of the GaN band gap 

1.5.3 The optical constants of GaN 

1.5.3. 1 The refractive index of GaN 

1.5. 3. 2 The absorption coefficient of GaN 

1.5.3 .3 Emission properties of GaN 7 

1.6. Properties of doped GaN 

1.6.1 Mg doping 

1.6.1 Cd Doping 

1.6.3 Be Dopig 

1.6.4 Hg Dopig 

1.6.5 Group IV doping 

1.6.6 Li doping 

1.7 Growth techniques 

1.7.1 Molecular beam epitaxy (MBE) 



1 

ii 

iii 
vi 

1 

1 

2 

Erreur ! Signet non defini. 

2 

Erreur ! Signet non defini.2 

2 

3 

3 

3 

4 

4 

5 

6 
6 

7 

8 
9 
9 

10 
10 
10 
10 
10 
11 



SUMMARY 



1.7.2 Halide vapor phase epitaxy (HVPE) 




11 


1.7.3 Metal-organic chemical vapor deposition (MOCVD) 




12 


Chapter II: The Physics of solar cells 




14 


II. 1 Introduction 




14 


II. 2 Solar energy 




14 


II.2.1 The Sun 




14 


II. 2. 1.1 The structure of the Sun 




15 


II.2.1.2Light 




15 


II.2.1.3Waves 




16 


II.2.1.4Absorption of light 


16 




II.2.2 Sunlight reaching earth 




17 


II. 2. 3 The impact of the movement of the earth on the solar radiation 


19 


II.2.4 Ground measurements of solar radiation 




19 


II.2.4.1 Pyrheliometer 




19 


II. 2. 4. 2 Pyranometer 




20 


II.3 p-n Junctions 




21 


II.3.1 Structure and principle of operation 




21 


II. 3. 2 Flat band diagram 




22 


II.3.3 Thermal equilibrium 




23 


II. 3. 4 Forward and reverse bias 




24 


II.4 Physical processes in p-n junction of solar cells 




25 


II.4.1 Equivalent circuit of the solar cell 




27 


II.4.2 Solar cell parameters 




28 


II.4.2.1Short circuit current (I sc ) 




29 


II.4.2.20pen circuit voltage (V oc ) 




29 


II.4.2.3FU1 factor (FF) 




29 


II.4.2.4Efficiency (q) 




29 


II.4.2.5Series resistance (R s ) 




30 


II.4.2.6Shunt resistance (R s h) 




30 


II.4.2.7Quantum efficiency (QF(X)) 




30 


II.4.3 Multijunction solar cell 




30 


II.4.3.1Bandgap Selection 




31 



IV 



SUMMARY 



II. 5 Semiconductors 




32 


II.5.1 Band Theory of Solids 




33 


II. 5. 2 Recombination and generation 




33 


II.5.2.1Band-to-band generation 




34 


II.5.2.2Band-to-band recombination 




34 


n.5.2.3 R-G center recombination orindirect thermal recombination 


34 


II. 5. 3 Equilibrium carrier concentrations 




35 


Chapter III: The SILVACO ATLAS Simulator 


38 


in.l Introduction 




38 


III. 2 Silvaco Atlas Simulation Software 




38 


III.2.1 Atlas Inputs and Outputs 




38 


III.2.2 Using "Silvaco- Atlas" software 




39 


HI.2.3 Running Atlas inside Deckbuild 




40 


III. 2. 4 Input file structure 




40 


III.3 Structure Specification 




41 


III.3.1 Mesh 




41 


III. 3. 2 Region 




42 


III.3.3 Electrodes 




43 


III. 3. 4 Doping 




44 


III.3.5 Materials Model Specification 




44 


III. 3. 5.1 Material 




44 


III.3.5.2Models 




45 


III.3.5.3Contact 




45 


III. 3. 5. 4 Interface 




45 


ID. 3. 6 Numerical Method Selection 




46 


m.3.7 Solution Specification 




46 


III.3.7.1 Log 




46 


III.3.7. 2 Solve 




47 


m.3.7.3 Load and Save 




47 


III. 3. 8 Results Analysis 




47 


Chapter IV: Results and discussions 




48 


IV. 1 Introduction 




48 


IV. 2 Model description 


4Erreur ! 


Signet non defini. 



V 



SUMMARY 



IV. 3 Simulation parameters 4Erreur ! Signet non defini. 

IV .4 Results and discussion 51 

IV. 4.1 Results for proportion of indium x=0 (GaN) 52 

IV. 4 . 2 Results for proportion of indium x=0.20 (Ino.2Gao.8N) 53 

IV. 4 . 3 Results for proportion of indium x=0.39 (Ino.39Gao.6iN) 55 

IV. 4. 4 Results for proportion of indium x=0.52 (Ino.52Gao.48N) 56 

IV. 4 . 5 Results for proportion of indium x=0.57 (Ino.57Gao.43N) 57 

IV. 4 . 6 Results for proportion of indium x=0.65 (Ino.65Gao.35N) 59 

IV. 4 . 7 Results for proportion of indium x=0.68 (Ino.68Gao.32N) 60 

IV. 4 . 8 Results for proportion of indium x=0.78 (Ino.78Gao.22N) 62 

IV. 4 . 9 Results for proportion of indium x=0.84 (Ino.84Gao.i6N) 63 

IV. 4 . 10 Results for proportion of indium x=0.89 (Ino.89Gao.nN) 65 

IV. 4 . 1 1 Results for proportion of indium x=[0-0.89] 66 

IV. 5 Summary of indium content in the InGaN (x) absorber layer effect 67 

IV . 5 .lThe effect on the short circuit current 67 

IV . 5 .2 The effect on the open-circuit voltage V oc 67 

IV . 5 .3 The effect on the fill factor (FF) 68 

IV . 5 .4 The effect on the efficiency 69 

IV. 6 Comparison between best results of p-n and n-p at x=0.84, b^Ga^N 69 
Chapter V: General Conclusion 71 

General Conclusion 71 

REFERENCES 72 



VI 



General Introduction 



GENERAL INTRODUCTIN 

Nitride-based materials such as In x Gai_ x N have become important in fabricating 
photovoltaic devices due to their band gap energies lying between 0.7 and 3.4 eV [1- 
2]. It can absorb the full solar spectrum by a single material of In x Gai_ x N with 
different indium contents. Many theoretical calculations on the performance 
evaluation of InGaN-based solar cells have been performed [1-3]. They demonstrated 
that a very high solar efficiency could be achieved. The direct bandgap properties and 
high absorption coefficients make it have great potential in photovoltaic application. 
Further, the GaN-based materials have a high resistant to high energy irradiation and 
temperature variations [5, 6]. Thus, it is suitable for space applications or high 
concentrator solar cell systems. However, a good-quality InGaN film with high 
indium composition cannot be obtained due to the low miscibility of InN in GaN. 

The goal of our work isstudy the effect of the proportion ofindium of the GaN 
absorber layer on the density of current in p- InGaN/ n-InGaN solar cell. For this we 
have used numerical simulation to study this effect. The numerical simulation is based 
on the SILVACO TCAD software. 

This memory is organized in a general introduction, four chapters and a conclusion. In 
the introduction the importance of the work is briefly outlined. In the first chapter, we 
introducedthe GaNmaterial In particular its chemical, crystalline andoptical 
propertiesare detailed. In the second chapter, we present the solar radiation. The 
operation ofthe p-n junction and the solar cell are also discussed. The third chapter is 
dedicated to give a detailed description on the software (SILVACO TCAD) used for 
the stimulation of our device. Finally, the fourth chapter presents and discusses the 
numerical results obtained of the simulation of a p-In x Gai_ x N/n- In x Gai_ x N solar cell. 
The work finished with a conclusion which summarizes the realized work. 
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Chapter I 

Properties of Gallium Nitride 



1.1 Introduction 

As all other Ill-nitrides Gallium Nitride (GaN) has partially covalent and partially ionic 
bonds. The ionic bond is relatively higher that explains the significant hardness and chemical 
stability of GaN. Other superior material properties such as high electron mobility, highthermal 
conductivity, and high critical field combined with direct wide band gap make GaN attractive for 
many electronic and optoelectronicapplications. 

I.2Crystallographic Properties 

GaN crystallizes in two different forms at room temperature: the hexagonal or wiirtzite structure 
and cubic structure or zincblende (sphalerite) (ZnS). The morphological forms of GaN crystals 
of the hexagonal structure is shown in Figure 1.1. 
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Figurel.l: Two morphological forms of GaN crystals: (a) Hexagonal platelets et (b) Hexagonal 

needels. Grid: 1mm 



At high pressure, it adopts a NaCl-type structure. These different structures are shown in Figure 

1 . 2 . 
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Figurel.2: The crystallinestructuresofGaN: (a) wiirtzite, (b) Zineblendeand(c) NaCl. 
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I.2.1Thehexagonalorwtirtzitestructure 

Thisstructure is also known as a-GaN. It isaverythermodynamicallystable structure, constituted 
bytwocompacthexagonalsubnets, one containing thegalliumatoms(Ga) and the other 
nitrogenatoms(N), andinterpenetratingoffset alongthe axiscof5/8thof the unit cell. Thisstructure 
is characterized bythe lattice parametersa and c. In reciprocal space, the firstBrillouin zonehas a 
hexagonalsection. 

I.2.2Thestructureofzincblende typecubic 

Thisstructure is also known asp-GaNThermodynamicallymetastable, it consists 
oftwosubnetworks:face centeredcubicgalliumatomsand nitrogenoffset by onefourth ofthe 
diagonal of themesh.The reciprocalofthisnetworkstructureisaface-truncated octahedron. 

I.2.3The structure NaCltype 

This is thehigh-pressurewiirtzitestructurevariety,thetransitionbetween the twostructuresoccurs 
around50GPa. It is composed oftwo subface centeredcubic latticesof atomsof gallium 
andnitrogen, shifted by a/2. 

1.3 Chemical properties of GaN 

GaN wasfirstsynthesizedbyJohnsonet al[7] inl928 who noted that "Gallic nitride is an 
exceedingly stable compound". It is this chemical stability at elevated temperatures combined 
with its wide band gapthat has made GaN an attractive material for device operation in high 
temperature and caustic environments. GaN is also an excellent candidate for protective coatings 
due to its hardness. However, the majority of GaN research is interested in semiconductor device 
applications. While the thermal stability of GaN should allow high temperature processing steps 
to be performed, the chemical stability of GaN presents a technological challenge. There have 
been many reports describing the resistance of GaN films to the conventional wet etching 
techniques used in semiconductor processing. 

1.3.1 The surface chemistry of GaN 

Several groups have studied the surface chemistry of GaN usingthe techniques of Auger 
electronspectroscopy (AES)[8], x-ray photoemissionspectroscopy (XPS)[9], andelectron energy- 
loss spectroscopy (EELS) [8]. There were also a number ofstudies examining 
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thethermalstabilityanddissociationof GaN.However, the materialstudied wasgenerallynot 
wellcharacterized and the results from different laboratories varied widely. 

1.3.2 The stability of GaN at hightemperature 

Studiedthe stability of GaN at hightemperature, significant weightlosseswerereported 
attemperatures aslowas750 ° C byFurtado and Jacob[10] 

whilenosignificantlosseswereobservedatlOOO°C by Morimoto[l 1], 

ObservedthatGaNwas less stable in a HC1 or FFatmosphere while the best results were 
obtainedinN2[10]. Somecontroversyexists overthereaction whichdominatesthe decompositionof 
GaN. MunirandSearcy[12]measuredthe apparent vapor pressure of GaN and calculatedthe heat 
of sublimation to be 72.4+0. 5kcal/mol. 

1.4 Electrical properties of GaN 

Control ofthe electrical properties of GaN remainsthe foremost obstacle hindering device efforts. 
Unintentionally doped GaN has in all cases been observed to be n type with 
thebestsamplesstillshowinganelectronconcentration n ~ 4 X10 16 /cm 3 [8]. No impurity has been 
shownto be present in sufficient quantityto account for the carriers, so 

researchershaveattributedthebackgroundtonativedefectswhichare 

widelythoughttobenitrogenvacancies. Allefforts, untilrecently,toobtainp-typedopinghave resulted 
in heavily compensated, highly resistive films. The development of a reliable p-type doping 
technique for GaN remains aprimary challenge forresearchers. 

1.4.1 GaN mobilities 

Nakamura [13] hasreportedacceptable GaN mobilitiesin 1991. Room- 
temperatureandliquidnitrogentemperature mobilitiesof p n ~ 600cm /V sandp n ~ 1500 cm ~/V 
satcarrier concentrationsofn = 4 X 10 16 /cm 3 andn=8X10 15 /cm3 weremeasured.Ilegemsand 
Montgomery [ 1 4]performedtemperaturedependentHall measurementsfrom 1 0- 

900KonGaNhavingaroom-temperaturebackgroundelectronconcentrationn=3 X 1 0 /cm 

andaroom-temperaturemobilityofpn =~ 300 cnr/Vs. Their data are shown in FigureI.3.In their 
best layer, the mobility peaked at p n =~ 800 cm 2 /Vsat 150 Kand decreased atlowertemperatures 
whichthe authors attributedto the onset of impurity band conduction. Inthe 300-900 Krange,the 
mobilitytended to vary with T 2 ,and at 900 K, a mobility ofp n =~ 25 cm 2 /Vs was measured. 
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Figurel.3: Temperature dependence of the electron mobilityinunintentionally doped GaN 
havingroom-temperature background electronconcentrations(U43Al) 
n=2xI0 17 /cm , and(U51Al)n=4xl0 17 /cm J [14]. 

I.50ptical properties of GaN 

GaNisprimarilyof interestforitspotentialasablue and UVlight emitter.For thatreason, much of the 
effort-devotedto GaN has been directedtowards determiningitsoptical properties. A fairly 
complete account ofthe optical properties of wurtzite GaNispossible. 

1.5 .1 Band gap energy ofGaN 

Maruska and Tietjen[15] were the firstto accurately measurethe GaN directband gap energyto be 
3.39 eV(FigureI.4). Soon after, Pankove et al[ 16] reported a lowtemperature (1.6 K) GaN PL 
spectrum. Strong near gap emission was observed at 3.477 eVand a weaker peak was noted at 
3.37 eVwhichsplitintoa~10 meVspaced doublet whenthe temperaturewasincreasedto35K.The 
separation ofthe two bands wastaken torepresentalower boundof thedonor bindingenergy. Using 
a value of () =9.8 obtained fromthe GaN phonon energies, andadonorbindingenergyof 42+1 
meV,an electroneffective mass m* = 0.25m o wascalculated.An acceptor binding energy of -200 
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meV was also estimated. At 300 Kthe energy of the highestenergytransition had decreasedto 3.39 
eV whichisin excellent agreement with the observations of Maruska and Tietjen[16]. 




Figurel.4: Room-temperature absorption spectrum of GaN [11]. 

1.5.2 Thetemperatureand pressuredependenceoftheGaNbandgap 

SeveralgroupshavemeasuredthetemperaturedependenceoftheGaNbandgap.Pankoveetal[16]estima 
ted abandgaptemperaturecoefficientofdE g /(dT) = — 6.0xl0 4 e V IK inthe linear region above 
180 °C fromthetemperature dependence of the principle emission peak.Monemar[17] 
determinedthefundamental bandgaptobe3.503 ± 0.0005eVatl.6Kandfitedthe temperature 
dependenceof thebandgaptotheempirical relation: 

Eg = 3.503+(5.08 x 10' 4 T 2 ) / (T-996) eV (1.1) 

Opticaltransmissionmeasurementswere also used tostudy its evolutionas a functionof pressure 
[18]: 

E g = Eg(P=0) + 0,47.P + 0,0018.P 2 eV(I.2) 

Camphausen and Connell [19] studied the pressure and temperaturedependenceof the 
GaNabsorptionedge. Under hydrostaticpressureupto 10 kbar,theGaNbandgap 

energyshiftedbydE g /(dP)=4.2 ± 0.4 meV /kbarat roomtemperature. The temperature dependence 
was measuredtobe dE/(dT) = - 6.7x10 4 eV/K. 
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1.5.3 The optical constants of GaN 

In these studies. Although there have been many reports on the optical properties and constants 
of GaN, the current state of the available data is seemed to be less than adequate. In particular, 
the optical constants. 

1.5.3. 1 Therefractiveindexof GaN 

There have been several measurements oftheoptical constants of GaN. Ejder [20] 
measuredtherefractiveindexof GaNusingtransmissionandabsorption 

measurementsandobtainedvaluesof n(3.38eV) =2.67 andn(l.OeV) = 2.33. 

In the particularcase where the materialis isotropic, thecomplex refractive index 
is writtenas follows: 



wherenisthe real part ofthe refractive indexandthe extinction coefficientK. 
Thecomplex indexis related to thedielectric constantbythe relationship: 

r = n 2 = n 2 — K 2 + 2inK(I.4) 



such as: r = 1 + iz(L5) 

It comes = n 2 — K 2 and 2 = 2nK(I.6) 

In the range oftransparencyaway from thegap energyof the semiconductorband, 
thevariationofthereal indexofrefraction as a functionofthe wavelengthis given by 
theSellmeierlaw [2 1 ] : 



N = n + iK 



(1.3) 




a, b et c are the Sellmeiercoefficients. This is shown in Figure 1.5. 
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D 




Wavelength, pm 

Figurel.5: Refractive indexofGaNbetween0.35andl0 jum wavelength. 



I.5.3.2 The absorption coefficient of GaN 

The absorption coefficienthasalsobeen the subjectofseveral studies. For materialstodirect gap, the 
value of the square ofthe coefficient islinearwith the energy ofphotonsabove theenergygap, that 
is: 



a 2 = (a 0 ) 2 (E-E g )(I.8) 

This coefficientis related to theextinctioncoefficientby the following equation: 

47TV 



a — 



K( 1 . 9 ) 



c and v are respectively the speed and frequency of the light. 



1.5.3 .3Emission properties of GaN 

Stimulated emission has been achievedin GaN by optical pumping in a number 
oflaboratories. Dingle et al[22] werethefirsttoreportlasinginGaN.Theirsingle crystal GaN needles 
emitted at aphoton energy of ~ 3.45 e Vat2K.Catalanoetal[13] wasabletoobserve80K 
stimulatedemissioninGaNat3.435eVwithanoptical gain of ~ 10 cm" .Amano et al[23] 

observedstimulated emissionatroomtemperatureforthefirsttime.Khanet al.[24] 

havesincereportedroom-temperaturestimulated emissionin GaN froma vertical(c-oriented)cavity. 
A summary of the structural, optical, properties of WiirtziteGaN is presented in Table 1.1. 
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Table 1. 1: Properties of wtirtziteGaN 



GaN 


Parameter value 


References 


Band gap energy 


Eg(300 K) =3.39 eV 


[16] 




Eg(1.6 K) =3.50 eV 


[18] 


Band gaptemperature 
coefficient(T> 1 80 K) 


AE8 /(dT) = ' 6 - 0X 10 ~ 3eV/k 


[17] 


Band gap pressure 
coefficient^ = 300 K) 


AEg /(dpy 4 - 2 x 10 ' 3eV/kbar 


[26] 


Lattice constant(T=300K) 


a =3.18910" 1() m 
c=5.18510' 10 m 


[16] 


Coefficient of thermal 


Aa / a =5.59 X10' 6 K _1 


[16] 


expansion(T = 300 K) 


Ac / a = 3.17X10' 6 K'‘ 




Thermal conductivity 


K= 1.3 W/cmK 


[27] 


Index ofrefraction 


n(l eV) =2.33 


[21] 




n(3.38 eV)=2.67 




Electron effective mass 


me=0.20 +0.02m o 


[28] 



1.6. Properties of doped GaN 

There has been much effort aimed at doping p type GaN byintroducing groupll andIV elements. 
Many potential dopants have been incorporated into GaN and their optical properties, shown 
schematically 

inFigureI.6.SomeimpuritieshavebeenobservedtoeffecivelycompensateelectronsinGaNleadingtohi 

ghlyresistivematerial. 
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Ga Site 



conduction Band 



30 met 



N Site 

" 

Native Defect Leve 



750mev Li 

700 mev Bt 

550 mev , d 

410 mev H „ 

250 mev Mg 



860 mev 



C 



225 mev 



Si 



Valence Band 

FigureI.6:Estimatedionisationenergies of somesubstitutionalacceptorsin GaN. Also shown isthe 
shallow donor level, generallythought to bethe nitrogen acancy, whichisresponsibleforthe n-type 
background carrier concentration observedin GaN[17]. 

1.6.1 Mg doping 

Akasakiand co-workers [28]developedatechniquewhichconverts Mg-compensated GaNinto p- 
type materialasaresultofelectronbeamirradiation.Nakamuraet al[29] havesucceededingrowingp- 

1 C O 

typeGaNbymetalorganicCVD (MOCVD)whichhadaholeconcentrationp=2 xlO /cm .Optical 
measurementshaveshownthatthe acceptor levelsin GaN are severalhundred milli-electron- 
voltsabovethe valence band edge. In order to grow GaN whichhas significant p-type conductivity 
at roomtemperature, existing Mg-doping methods mustbeimprovedor 

acceptorcomplexeswhichlieclosertothevalenceband edge must be sought. 

I.6.2Cd Doping 

Therehavebeenseveralreportswhichdescribeattemptstodope GaN withCd.Bergmanet al[30], were 
ableto dope GaN withasmall amount of Cd,but a 

sufficientamounttocompensatetheelectronbackground 

concentrationcouldnotbeincorporated.Ilegemsetal[31]reportedtheCd-related peak at 2.85 
eV.Bergman et al [30], measured alps decaytime at77Kfortheir 2.7eVemission. 
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I.6.3Be Dopig 

Twogroups reportedBe-dopedGaN.Bewas observedtoeffectivelycompensatethefreeelectronsand 
highresistivity GaN wasobtained. Optical measurements detected a broad Be-related emission at 
an energy of ~ 2.2 eV. 

I.6.4Hg Dopig 

Two groupshaveinvestigated Hg doping.No electricalmeasurementswerereported.Pankoveand 
Hutchby [32]observed2.43eVemissioninHg-dopedsamples. Ejder and Grimmeiss [33] observed 
2.9 eV emission and estimatedthe Hg-acceptorleveltobe410 meVabovethe valence band. 

I.6.5Group IV doping 

Group IVelementshave beeninvestigatedaspossible p-typedopantsinthehopesthattheywould 
occupy vacant mtrogen sites in GaN. Cunningham et al[34]observedanorder of 
magnitudereductioninthe freecarrierconcentrationwhenGedopingwasattempted. Maruska and 
Tietjen [15] reported occasional p-type conductionin Ge-doped GaN, but laterretracted this claim 
[13]. In thesamework,Siwasobservedtomarkedlyreducethe 

electronconcentrationinGaNtoaslittleasn=3 x 10 16 /cm 3 butp-typeconductivitywasnotachieved. 

I.6.6Li doping 

Anumberof workers haveinvestigatedLi doping. Pankove et al [35] obtained compensated GaN 
using Li doping. The Liacceptorluminescence wasobservedat2.23eV.Ejder and Grimmeiss [33] 
estimated the Li-acceptor levelto be 750 meVabovethe valence band. 

I.7Growth techniques 

When heated at atmospheric pressure GaN dissociates rather than melts, thus growth from liquid 
solution at reasonable conditions is impossible as it mentioned above. Therefore vapor phase 
deposition techniques are used for GaN synthesis; three major methods exist nowadays, MBE, 
HVPE and MOCVD, which will be described in more detail. 
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1.7.1 Molecular beam epitaxy (MBE) 

MBE is a non-equilibrium vapor phase epitaxy growth method where solid or gas source 
elements are heated in evaporators, known as Knudsen effusion cells, to produce beams of atoms 
impinging on the substrate, which is rotated and heated to appropriate temperatures. The effused 
atoms need to have mean free paths bigger than the size of the chamber, also to obtain films of 
sufficient purity low levels of residual gases must be sustained, therefore ultra- high vacuum (10 9 
Torr) is needed, which is the drawback of this method. The advantages are the lower, compared 
to other methods, growth temperatures of 600 - 800 °Cand possibility to grow very thin layers 
with fine control of composition. Precision of this growth technique is based on low deposition 
rates and in-situ characterization capabilities. For example, reflection high energy electron 
diffraction (RHEED) detector allows to control the thickness of the film up to one atomic 
monolayer. 



Vacuum 

pump 



mass spectrometer 

RHEED gun 



Figurel.7: Schematic drawing of MBE chamber for GaN growth. 

1.7.2 Halide vapor phase epitaxy (HVPE) 

A typical horizontal hot wall HVPE reactor is depicted on Figurel.8. The chamber is made of 
quartz and can be heated resistively or by RF induction. There are two temperature zones. The 
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metallic source zone is kept at the temperature of- 800-900 °C, which is lower than the growth 
zone temperature (1000-1 100°C). 

By using several temperature zones more flexibility in controlling and in adjustingGrowth 
conditions can be achieved. 



800-900 °C DC . lc 1000-1100 9 C 




Figurel.8: Schematic drawing of HYPE growth reactor. 



Hydrogen chloride is let into the quartz tube with melted Ga where the reaction producing 
gallium chloride occurs: 

Ga (liq) +HCl( g ) - GaCl (g) +^H2( g) (1. 11) 

The efficiency of this processis very high, above 95%. Inaseparate tube NH 3 , the group V source, 
is flowed. Subsequently ammonia and gallium chloride undergo the following reaction on the 
surface of the substrate: 

GaCl (g) +NH 3(g) -> GaN+HCl (g )+ H 2(g) (1.12) 

To prevent parasitic growth in the gas inlet ammonia and halide (GaCl) should not be mixed 
before they reach the substrate. Light molecule gases such as H 2 or mixtures of N 2 and 
H 2 areusedascarriergasestoobtainlaminar flow.The deposition rate is - 100- 300 pm/h and is very 
highcompared to MBE, which is why HVPE is commercially used for growing thick GaN layers. 
The process occurs at atmospheric pressure, at near equilibrium conditions and is mass transport 
limited by the flow of halide precursor. 

1.7.3 Metal-organic chemical vapor deposition (MOCVD) 

MOCVD is a technique similar to HVPE, only instead of inorganic species metal-organic 
precursors are flown in the reactor as sources of Ga. Specifically trimethyl-gallium ((CH 3 ) 3 Ga) is 
used as Ga, and ammonia-as nitrogen precursors. Organic species and ammonia are transported 
in separate pipes to the quartz chamber to prevent premature reacting. Lighter carrier gases (H 2 
or N 2 ) are let through organic compound liquids thus transporting the precursors to the reactor. 
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These compounds are thermally dissociated over the surface of ~ 1000°Chot substrate and react 
in a sequence of complex reactions forming GaN. In hot wall MOCVD the whole reactor is 
heated by RF coils, whereas in the cold-wall process only the susceptor, which is usually made 
of graphite, is heated. The substrate may be rotated to obtain uniform films. Before the growth is 
started nitridation of the substrate is usually done to reduce lattice mismatch due to the nitrides 
formedon the substrate surface (SiN on SiC or AIN on AFO 3 ). The usual practice is to start 
epitaxial growth on sapphire substrates with a low-temperature (~ 600°C) GaN buffer layer. Ease 
of scalability, high uniformity and quality of the films, and a relatively high growth rates 
(~lpm/h) are making this technique the method of choice for commercial purposes. 



Page 13 



Chapter II 



The Physics of solar cells 



Chapter II 

The Physics of solar cells 



II. 1 Introduction 

Semiconductor solar cells are fundamentally quite simple devices. Semiconductors have the 
capacity to absorb light and to deliver a portion of the energy of the absorbed photons to carriers 
of electrical current - electrons and holes. A semiconductor diode separates and collects the 
carriers and conducts the generated electrical current preferentially in a specific direction. Thus, 
a solar cell is simply a semiconductor diode that has been carefully designed and constructed to 
efficiently absorb and convert light energy from the sun into electrical energy. 

II.2 Solar energy 

The sun is responsible for nearly all of the energy available on earth. The exceptions are 
attributable to moon-tides, radioactive material, and the earth's residual internal heat. Everything 
else is a converted form of the sun's energy: Hydropower is made possible by evaporation- 
transpiration due to solar radiant heat, the winds are caused by the sun's uneven heating of the 
earth's atmosphere, fossil fuels are remnants of organic life previously nourished by the sun; and 
photovoltaic electricity is produced directly from sun-light by converting the energy in sunlight 
into free charged particles within certain kinds of materials. 

Solar energy can be used to produce electricity by converting solar radiation directly into an 
electric current through the so-called solar cells or photovoltaic cells. The electrical power 
produced can either be injected directly into the grid (grid connected) or used in-situ (off-grid). 

II.2.1 The Sun 

Schematic structure of the Sun is shown in Figure II.1. The Sun’s energy is the result of 
severalnuclear fusion reactions, the main being the process in which hydrogen (fourprotons) 
merges and forms helium. The mass of alpha particle (helium nucleuses less than the mass of 
four protons, the difference in mass is converted into energy according to Einstein’s formula 

E = me 2 (II.l) 
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II.2.1.1The structure of the Sun 

The Sun is a large sphere of intensely hot gases consisting, by mass, about 75 % of hydrogen, 23 
% of helium and others (2 %). This proportion changes slowly over time referring to the nuclear 
fusion in its core with temperatures of approximately 15-20 million K. Hydrogen atoms fuse 
there to form helium and this energy is then delivered as radiation (light and heat) into space. 
The Sun’s outer surface, namely photosphere, has an effective blackbody temperature of approx. 
6000 K. This mean, as viewed from the Earth, the radiation emitted from the Sun appears to be 
essentially equivalent to that emitted from a blackbody at 6000 K [37]. To understand the 
behaviour of the radiation from the Sun the characteristics of the blackbody should be discussed 
here. 
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Figure II.l:The structure of the Sun. 



II.2.1.2Light 

The Earth receives energy from the Sun by way of radiated energy. This radiated energy is 
carried by tiny particles having no mass called photons. Photons behave like waves and, as such, 
have a characteristic wavelength, frequency, and energy. It is the energy of these photons that is 
used in photovoltaic cells to excite electrons so that an electrical current can be produced. 



Page 15 



Chapter II 



The Physics of solar cells 



II.2. 1.3 Waves 

Figure II.2illustrates that the wavelength of a wave is the distance between two successive peaks, 
or maxima, of a wave. 



Wavelength 

4 — ► 




Figure II.2:Wavelength of wave. 

The frequencyof a wave is a measure of how many peaks are observed in one second froma 
stationary point along a wave’s path. The product of the wavelength and frequency gives the 
wave’s velocity: 

v = If (11.2) 

All radiative, or light, energy travels at the same velocity, commonly known as the speed of 
light, which is about 300,000,000 m / s (3 x 10 8m / s ) Because the velocity of light is a constant 
and also a product of frequency and wavelength, specific frequencies of light are associated with 
specific wavelengths. 

II.2.1.4Absorption of light 

Photons incident on the surface of a semiconductor will be either reflected from the top surface, 
will be absorbed in the material or, failing either of the above two processes, will be transmitted 
through the material. For photovoltaic devices, reflection and transmission are typically 
considered loss mechanisms as photons which are not absorbed do not generate power. If the 
photon is absorbed it has the possibility of exciting an electron from the valence band to the 
conduction band. A key factor in determining if a photon is absorbed or transmitted is the energy 
of the photon. Therefore, only if the photon has enough energy will the electron be excited into 
the conduction band from the valence band. Photons falling onto a semiconductor material can 
be divided into three groups based on their energy compared to that of the semiconductor band 
gap: 

E p h< Eg Photons with energy E p h less than the band gap energy Eointeract only weakly with the 
semiconductor, passing through it as if it were transparent. 

Eph = Eg have just enough energy to create an electron hole pair and are efficiently absorbed. 



Page 16 



Chapter II 



The Physics of solar cells 



E p h> Eg Photons with energy much greater than the band gap are strongly absorbed. However, 
for photovoltaic applications, the photon energy greater than the bandgap is wasted as electrons 
quickly thermalize back down to the conduction band edges. 

II.2.2 Sunlight reaching earth 

Even though the sun ranks as a run-of-the-mill star, it releases a huge quantity of energy in terms 
of human capacity or need. Power output per second is 3.86 x 10" megawatts (MW). This 
energy fills the solar system, bathing the earth's atmosphere with a near constant supply of 1.37 
kilowatts per square meter (kW/m ).Not all of the direct sunlight incident on earth's atmosphere 
arrives at the earth's surface. The atmosphere attenuates many parts of the spectrum (Figure II.3). 
For example, X-rays are almost totally absorbed before reaching the ground. A good percentage 
of ultraviolet radiation is also filtered out by the atmosphere. Some radiation is reflected back 
into space. Some is randomly scattered by the atmosphere, which makes the sky look blue. 




Figure II.3:Spectral power density of black-body radiation at 6000 K, AMO and AM1.5 

spectra. 

It is valuable to relate the amount of sunlight at the earth's surface to the quantity, or air mass 
(AM), of atmosphere through which the light must pass. Radiation arriving at the surface of the 
earth is measured against that reaching the fringes of the atmosphere, where there is no air, and 
the air mass is zero (AMO). The light of the high-noon sun (and under further specified 
conditions) passes through an air mass of one (AMI). The intensity of the sunlight reaching the 
ground weakens for sun angles approaching the horizon since the rays have more atmosphere, or 
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air mass, to penetrate. The atmosphere is a powerful absorber and can cut the sun's energy 
reaching the earth by 50% and more. 




Figure II.4:The Path length (in unit of air mass) changes with the zenith angle. 

The Earth revolves in orbit around the Sun, that is approximately circular circumference 
deviation does not exceed 1.7%. The average distance between the Sun and the Earth (Figure 
H.5), called the astronomical unit is equal to 1,495 1011 m. The Sun, as viewed from the Earth’s 
surface, represents a disk with an angular diameter 32' 



Sun Earth 
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II.2.3 The impact ofthe movement ofthe earthon thesolar radiation 

The Earth rotates in an elliptical orbit around the Sun during 1 year with the Sun in one of the 
focal points. Also, the Earth rotates around its axis, making one revolution in 24 h. The Earth’s 
axis of rotation has a fixed direction in space and is inclined at angle 5o= 23.44° towards the 
normal line in the orbit plane (Figure II. 6). 




Figure II.6:Earth’s orbit and the declination angle 5. 



The angle between the Sun direction and the equatorial plane is called the declination angle 5 
and varies throughout the year from +23.44° when the summer solstice (June 20/21) to 
— 23.44°at the winter solstice (December 21/22). On March 21st, and respectively, September 
22/23 the declination 8= 0. 

II.2.4 Ground measurements of solar radiation 

Radiometry is the science of electromagnetic radiation measurement. The generic device is 
named radiometer, for instance, the pyrheliometer that measures the direct beam irradiance and 
the pyranometerthat measures the horizontal beam and diffuse irradiances. 

II.2.4.1 Pyrheliometer 

The pyrheliometer is a broadband instrument that measures the direct beam component G n of 
solar radiation. Consequently, the instrument should be permanently pointed toward the Sun. A 
two-axis Sun tracking mechanism is most often used for this purpose. The detector is a multi- 
junction thermopile placed at the bottom of a collimating tube (Figure II.7a) provided with a 
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quartz window to protect the instrument. The detector is coated with optical black paint (acting 
as a full absorber for solar energy in the wavelengths range 0.280-3 pm). Its temperature is 
compensated to minimize sensitivity of ambient temperature fluctuations. The pyrheliometer 
aperture angle is 5°. Consequently, radiation is received from the sun and a limited circumsolar 
region, but all diffuse radiation from the rest of the sky is excluded. A readout device is used to 
give the instant value of the direct beam irradiance. Its scale is adapted to the sensitivity of the 
particular instrument in order to display the value in SI units, Wm"".For illustration, a picture of a 
Hukseflux DR01 First Class pyrheliometer (Hukseflux 2012) is presented in Figure II.7b. 




Figure II.7:a Schematic of a pyrheliometer. b Photo of Hukseflux DR01 first class 
pyrheliometer (Hukseflux 2012). (Public license on Wikimedia commons). 



II.2.4.2Pyranometer 

Pyranometers are broadband instruments that measure global solar irradiance incoming from a 
In solid angle on a planar surface. A typical pyranometer is schematically represented in Figure 
II.8a. It consists of a white disk for limiting the acceptance angle to 180° and two concentric 
hemispherical transparent covers made of glass. The two domes shield the sensor from thermal 
convection, protect it against weather threat (rain, wind, and dust) and limit the spectral 
sensitivity of the instrument in the wavelength range [0.29-2.8 pm]. A cartridge of silica gel 
inside the dome absorbs water vapor. 

A pyranometer can be also used to measure the diffuse solar irradiance Gd , provided that the 
contribution of the direct beam component is eliminated. For this, a small shading disk can be 
mounted on an automated solar tracker to ensure that the pyranometer is continuously shaded. 
Alternatively, a shadow ring may prevent the direct component Gb from reaching the sensor 
whole day long (see Figure II.8b). Because the daily maximum Sun elevation angle changes day 
by day, it is necessary to change periodically (days lag) the height of the shadow ring. 
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Figure II.8:(a) Schematic of a pyranometer. (b) First class pyranometer LPPYRA 12 
(DeltaOHM 2012) equipped with shadow ring, mounted on the solar platform of the West 
University of Timisoara, Romania (SRMS 2012). 



II.3 p-n Junctions 

P-n junctions consist of two semiconductor regions of opposite type. Such junctions show a 
pronounced rectifying behavior. They are also called p-n diodes in analogy with vacuum tubes. 
The p-n junction is a versatile element, which can be used as a rectifier, as an isolation structure 
and as a voltage-dependent capacitor. In addition, they can be used as solar cells, photodiodes, 
light emitting diodes and even laser diodes. 

II.3.1 Structure and principle of operation 

A p-n junction consists of two semiconductor regions with opposite doping type as shown in 
Figure II. 9. The region on the left is p-typc with an acceptor density N a , while the region on the 
right is /i-type with a donor density N d . The dopants are assumed to be shallow, so that the 
electron (hole) density in the / 2 -typc (p-type) region is approximately equal to the donor 
(acceptor) density. 
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x d 

p-type < n-type 




Figure II.9:Cross-section of a p-n junction. 

We will assume, unless stated otherwise, that the doped regions are uniformly doped and that the 
transition between the two regions is abrupt. We will refer to this structure as an abrupt p-n 
junction. 

Frequently we will deal with p-n junctions in which one side is distinctly higher-doped than the 
other. We will find that in such a case only the low-doped region needs to be considered, since it 
primarily determines the device characteristics. We will refer to such a structure as a one-sided 
abrupt p-n junction. 

The junction is biased with a voltage V a as shown in Figure II. 9. We will call the junction 
forward-biased if a positive voltage is applied to the p-doped region and reversed-biased if a 
negative voltage is applied to the p-doped region. The contact to the p-type region is also called 
the anode, while the contact to the n-type region is called the cathode, in reference to the anions 
or positive carriers and cations or negative carriers. 

II.3.2 Flatband diagram 

The principle of operation will be explained using a thought experiment, an experiment, which is 
in principle possible but not necessarily executable in practice. We imagine that one can bring 
both semiconductor regions together, aligning both the conduction and valence band energies of 
each region. This yields the so-called flatband diagram shown in Figure II. 10. 
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type and p- type regions. 

Note that this does not automatically align the Fermi energies, Ep, n and E|. p . Also, note that this 
flatband diagram is not an equilibrium diagram since both electrons and holes can lower their 
energy by crossing the junction. A motion of electrons and holes is therefore expected before 
thermal equilibrium is obtained. The diagram shown in Figure II. 10 (b) is called a flatband 
diagram. This name refers to the horizontal band edges. It also implies that there is no field and 
no net charge in the semiconductor. 

II.3.3Thermal equilibrium 

To reach thermal equilibrium, electrons/holes close to the metallurgical junction diffuse across 
the junction into the p-type/n-type region where hardly any electrons/holes are present. This 
process leaves the ionized donors (acceptors) behind, creating a region around the junction, 
which is depleted of mobile carriers. We call this region the depletion region, extending from x = 
-x p to x = x n . The charge due to the ionized donors and acceptors causes an electric field, which 
in turn causes a drift of carriers in the opposite direction. The diffusion of carriers continues until 
the drift current balances the diffusion current, thereby reaching thermal equilibrium as indicated 
by a constant Fermi energy. This situation is shown in Figure II. 11: 
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Figure 11.11: Band Diagram of p-n Junction in Thermal Equilibrium [38]. 

II.3.4 Forward and reverse bias 

We now consider a p-n diode with an applied bias voltage, V a . A forward bias corresponds to 
applying a positive voltage to the anode (the p-type region) relative to the cathode (the n-type 
region). A reverse bias corresponds to a negative voltage applied to the cathode. Both bias modes 
are illustrated with Figure 11.12. The applied voltage is proportional to the difference between the 
Fermi energy in the n-type and p-type quasi-neutral regions. 

As a negative voltage is applied, the potential across the semiconductor increases and so does the 
depletion layer width. As a positive voltage is applied, the potential across the semiconductor 
decreases and with it the depletion layer width. The total potential across the semiconductor 
equals the built-in potential minus the applied voltage, or 

- v a (ii.3) 
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Figure II.12:Energy band diagram of a p-n junction under reverse and forward bias. 



II.4Physical processes in p-n junction of solar cells 

During the illumination of the cell (junction) a number of physical processes schematically 
illustrated in Figure II.13will occur: 

1) Absorption of photons with energy hv> E g 

2) Generation of charge carrier pairs-electron and hole 

3) Migration of charge carriers (diffusion). 

4) Separation of charge carriers induced by built-in electric field. 

The charge separation results in emerging an electric voltage on the electrodes of the solar cell. 
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Figure II. 13: Electron processes in the region of p-n junction. A- absorption of photons, M - 
migration of charge carriers, S - separation of electron and holes, R - recombination. 



After the generation of the charge carriers (before the separation phase) the reverse and a non 
advantageous process - a recombination of electrons and holes could appear. Usually it is 
performed via defect centers. The recombination can happen in a volume of semiconductors 
constituting the PV structure as well as in the range of the surface. The surface recombination 
can be intensive and measures of reduction of these phenomena must be undertaken. 

A separation effect is the final one and is the most influencing the efficiency of the PV structure 
and the solar cell. The effect will be the most efficient if absorption, generation of charge carriers 
(or quantum yield) and migration of carriers to the interface range are high and advantageous. A 
high level of scientific knowledge and availability of advanced technologies are needed in order 
to fulfill this requirement. 

Absorber illuminated by the solar light absorbs photon when their energy is higher than the 
energy gap of the absorber. Passing distance x in a semiconductor the photon flux is absorbed 
according to the dependence 

NphO) = Np h (0) (1 - r)exp[-a(X)x] (II.4) 

Were r and a are coefficients of reflection and absorption, respectively. 

Figure 11.14 shows a plot of a vs (hv) for several semiconductors used in fabricating solar cells. 
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Figure 11.14: Absorption edges of the photovoltaic semiconductor. 



II.4.1Equivalent circuit of the solar cell 

Figure 11.13 represents the equivalent circuit of a solar cell. The current through a p-n junction 
diode in dark is given by 

I = I sat (e^-l)(II.5) 

Where I sat is reverse saturation current, n is diode quality factor, V is applied voltage, K is 
Boltzmann constant and T is absolute temperature. 



Rs i 




Figure 11.15: Equivalent circuit of ideal solar cell and practical solar cell. 
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The total current flowing through the circuit is given by 



I=-I 



gen 




+ V ~ 1R s 
Rsh 



(II.6) 



Where R s series resistance, R sh shunt resistance and l gen is light generated current. 



II.4.2Solar cell parameters 

The current-voltage characteristics (I — V) for a solar cell under dark and light condition are 
shown in figure 11.16. The (I — V) curve is obtained by plotting the output current as a function 
of the voltage. 




The IV curve of a solar cell is the superposition of the IV curve of the solar cell diode in the dark 
with the light-generated current[39]. The light has the effect of shifting the IV curve down into 
the fourth quadrant where power can be extracted from the diode. Illuminating a cell adds to the 
normal "dark" currents in the diode so that the diode law becomes: 



I = Ic 



ex p (is) - A ~ *ph 



(II. 7) 



Wherelph light generated current and I s Saturation current, n ideality factor. 
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II.4.2.1Short circuit current (I sc ) 

The Short circuit current is found by setting V=0 in diode equation. This results in the following 
expression for I sc 

I sc = - Igen (II. 8) 

Where I gen is light generated current. 

The short circuit current is normlly as a current density, I sc (mA/cm ) which is given by 

Jsc = l f (II. 9) 

were, A is the area of the solar cell. 

II.4.2.20pen circuit voltage (V oc ) 

The open circuit voltage is given by 

Voc = A— Inf— + l) (11.10) 

9 V Isat ' 



II.4.2.3Fill factor (FF) 

Another characteristic property of solar cells of interest is the fill factor, FF which is defined as 



FF = 



IscWc 



I in V|n 

IscWc 



( 11 . 11 ) 



Were I m , V m are the current and voltage at maximum power, P n 



II.4.2.4Efficiency (q) 

The efficiency of a solar cell is deefined as the ratio of light generated power from the cell to the 
incident power. It is a measure of light energy successfully converted to an electrical energy. 
Efficiency q (%) is given by the relation. 



q = 



r out 



q = 



r in 

FFV oc J sc 



q = FF V oc J s 



( 11 . 12 ) 

(11.13) 

(11.14) 



2 

At 25°C, the standar reference spectrum irradiance is 100 mW/ cm". So, Pi n (Photon energy 
input). 
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II.4.2.5Series resistance (R s ) 

The current generated by photons is represented by an independent source. The two resistors are 
shown to model two losses in a solar cell. R s loss is primarily due to the contact resistance of the 
front and back contacts of the solar cell. Other component that contributes to series reistance is 
resistance in the semiconductor material. Typical values of R s are [1—5] Q cm". High R s values 
lead to a decrease in the solar cell which current in turn affects the efficiency. 

II.4.2.6Shunt resistance (R S h) 

The Shunt resistance is used to model leakage currents. The shunt resistance mainly arises due 
leakage current created within a solar cell. Some of the reasons for the occurrence of these 
leakage paths are due to recombination current and surface defects in the semiconductor. These 
leakage paths are not uniformly distributed and vary from one cell to another. The combined 
effect of all the leakage paths the term called shunt resistance R s h. 



II.4.2.7Quantum efficiency (QE(/,)) 

For ideal solar cell, every incident photon would generate an electron-hole pair and each of the 
photo carriers would make it to the depletion region where they could be separated and collected. 
Photons with energy less than the bandgap have insufficient energy to generate photo-carriers. 
Even if it has the sufficient energy, it need not contribute to the photocurrent. Quantum 
efficiency is defined as the ratio of the number of electons in the external circuit produed by an 
incident photon of given wavelengh. It is the measure of the effectiveness of a device to produce 
electronic charge from incident photons. Quantum efficiency describes response of the device to 
different wavelengths of light. 



QE(k) = - 



1 he l sc (X) 



q * p(V 



(11.15) 



Were q is the unit charge, X is the photon waavelength, h is Planck’s constant, c s the speed of 
light in air, I sc (k) is short circuit current and P is incident light power. 



II.4.3 Multijunction solar cell 

Multijunction cells use a combination of semiconductor materials to more efficiently capture a 
larger range of photon energies. They do so without sacrificing photovoltage or creating losses of 
heat to the degree of single-junction cells. 
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Figure II.17:A multijunction cell functions by absorbing sequentially lower energy light 

from the incident light. 

II.4.3.1Bandgap Selection 

Consider the spectrumofsunlight split and distributed over a variety of semiconductor materials 
by using a prism(see Figure 11.18). 



In that case, a different semiconductor material could be selected (see Figure II. 19a) that would 
best match each portion of the spectrum.Although such an arrangement is reasonable in theory, 
mechanically maintaining such a design would be problematic. 




White light 



Separated 

wavelengths 



Figure II.18:Refraction of light through a prism. 
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Figure II. 19: Alternate configurations for multij unction spectral splitting. 



A more viable solution is to stack, or grow, multiple layers of semiconductors with decreasing 
bandgaps (see Figure II. 19b). Top layers are designed to absorb higher-energy photons while 
transmitting lower-energy photons that are absorbed by lower layers of the cell. In theory, a 
number of these layers can be so arranged. This strategy is the characteristic feature of 
multijunction solar cells. 

Alloys of Group III and Group V elements, as well as other related compounds, lend themselves 
well to the design of multijunction cells. Indium Nitride (InN), Gallium Nitride (GaN), and the 
more commonly used GaAs are examples of such III-V materials. By carefully adjusting the 
compositions, a range of bandgap energies can be achieved. If a lower bandgap material were 
desired, less indium and more gallium would be used in the compound, up to the point at which 
no indium is included, and the resulting GaN would have a bandgap of about 3.4 eV. 



II.5 Semiconductors 

Understanding the fundamental theory behind the design of solar cells is very important for 
fabricating high efficiency solar cells. The main part of a solar sell is the p-n junction and the 
composition ofthe solar cell. These terms are explained in the following sections. 
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II.5.1Band Theory of Solids 

Atoms in a solid are closely spaced, and their outer orbitals influence one another. The result is 
that instead of required discrete energy levels, the outer energy levels become (smeared), 
forming bandsof allowable energy levels. At this point, the limitations ofthe Bohr model ofthe 
atomas a useful model are exceeded. At this point, electrons in these bands are not associated 
with any particular atom. To understand the electrical properties ofa material, it is helpful to 
consider the two bands commonly referred to as the valence bandand conduction band. The 
valence band encompasses the energy levels at which most of the electrons (and at absolute zero, 
allof the electrons) will be found. Electrons found in the higher-energy conduction band are 
considered free electrons. Free electrons are shared among all the atoms of the solid and are not 
associated with any one atom. Therefore, they are very mobile. Any applied electric field can 
create an electrical current with these electrons. 

For insulator and semiconductor materials, the valence and conduction bands are separated by a 
region called the bandgap. The bandgap is a region of forbidden energy levels, as defined by the 
quantumaspects of modern atomic theory. The difference between the lowest energy level of the 
conduction band and the highest energy level of the valence band is called the bandgap energy 
(Eg). 

At temperatures above absolute zero, electrons may possess enough thermal energy (at least Eg) 
to be excited across the forbidden energy gap into the conduction band. This occurrence is less 
likely if the bandgap energy is large, as in insulators. In conductors, the conduction band and the 
valence bands overlap, and conduction of electrons occurs, even at absolute zero. 




Valence band 

Insulator Semiconductor Conductor 



Figure II.20:Energy bandgaps. 

II.5.2 Recombination and generation 

When a semiconductor is disturbed from the equilibrium state the electron and hole 
concentrations change from their equilibrium values. A number of processes start within the 
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semiconductor after the perturbation in order to restore the equilibrium conditions. We call 
these processes recombination- generation (R-G) processes and they take care that a 
carrier excess or deficit is eliminated (if the perturbation is removed) or stabilized (if the 
perturbation is maintained). Generation is a process whereby electrons and holes are 
created. Recombination is a process whereby electrons and holes are annihilated or 
destroyed. 



II.5.2.1Band-to-band generation 

The generation of an electron-hole pair can occur when the energy to break a bond is 
delivered by a photon (photogeneration) or by heat (direct thermal generation). This process 
can be visualized in terms of the band diagram as the excitation of an electron from the 
valence band into the conduction band. 





Figure 11.21: photogeneration and direct thermal generation. 



II.5.2.2Band-to-band recombination 

This process involves the direct annihilation of a conduction-band electron with a valence-band 
hole and is often called direct thermal recombination. It is typically a radiative process, in which 
a photon is emitted that carriesan excess energy released in the recombination process. 
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Figure II.22:direct thermal recombination. 



II.5.2.3 R-G center recombination orindirect thermal recombination 

As the name suggests the annihilation or creation of electrons and holes does not occur directly 
but it is facilitated by a R-G center . The R-G centers can be special impurity atoms or 
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lattice defects. Their concentration is usually small compared to the acceptor or donor 
concentrations. The R-G centers introduce allowed energy levels (E T ) near the centre of the 
forbidden gap. An electron is trapped at the R-G center and consequently annihilated by a hole 
that is attracted by the trapped electron. Though this process seems to be less likely than the 
direct thermal recombination, it is at most operational conditions the dominant 
recombination-generation process in semiconductor.The process is typically non-radiative and 
the excess energy is dissipated into the lattice in form of heat. 



Heat 

(phonons) 

/ w v -* 



Figure II.23:R-G center (indirect thermal) recombination. 
II.5.3Equilibrium carrier concentrations 

When the semiconductor is in thermal equilibrium (i.e. at a uniform temperature with no external 
injection or generation of carriers), the Fermi function determines the ratio of filled states to 
available states at each energy and is given by 

f(E) = (E * Ep) (II.16) 

1+e 'kT 
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Where E F is the Fermi energy, k is Boltzmann’s constant, and T is the Kelvin temperature. As 
seen in Figure 11.21, the Fermi function is a strong function of temperature. At absolute zero, it is 
a step function and all the states below E F are filled with electrons and all those above E F are 
completely empty. As the temperature increases, thermal excitation will leave some states below 
E F empty, and the corresponding number of states above E F will be filled with the excited 
electrons. 

The conduction-band and valence-band effective densities of state (cm" ),Nc and 
N v , respectively, are given by 

3 / 

N c = 2 ( h2 " j (11.17) 

and 

* 3/ 

„ /2Trm„KT\ / 2 

N v = 2 ^ J j (11.18) 

When the Fermi energy, E F , is sufficiently far (>3kT) from either band edge, the carrier 
concentrations can be well approximated (to within 2%) as [40] : 

(e f -e c ) 

n 0 = N c e kt (11.19) 



Po= N v e kt (11.20) 

and the semiconductor is said to benondegenerate. In nondegenerate semiconductors, the product 
of the equilibrium electron and hole concentrations is independent of the location of the Fermi 
energy and is just: 

(Ey— E c ) -Eg 

n 0 p 0 = n? = N c N v e kt = N c N v e kt (11.21) 
where n; is the intrinsic carrier concentration. 

In an undoped (intrinsic) semiconductor in thermal equilibrium, the number of electrons in the 
conduction band and the number of holes in the valence band are equal; n 0 = p Q = ni, where ni is 
the intrinsic carrier concentration. The intrinsic carrier concentration can be computed from 
equation : 

(Ey— E c ) ~ e g 

n i = V N cN v e 2 kt = ^/NcNye 2xr(ll.22) 

The Fermi energy in an intrinsic semiconductor, E, = E F ,isgivenb: 

E i=^ + T ln 0 (IL23) 
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Under assumption that the semiconductor is uniformly doped and in quilibrium a simple 
relationship between the carrier and dopant concentrations can be stablished. We assume that at 
room temperature the dopant atoms are ionized. Inside a semiconductor the local charge 
density is given by: 

p = q (p+Nd - n - NJ) (11.24) 

where qis the elementary charge (q= 1.602 x 10" 19 C). Under equilibrium conditions, in the 
uniformly doped semiconductor the local charge is zero, which means that the semiconductor is 
everywhere charge-neutral. One can write: 

p + Ni - n - N2 = 0 (11.25) 

As previously discussed, the thermal energy available at room temperature is sufficient to 
ionise almost all the dopant atoms. Under this assumption,]^ ~ N D and N A ~ N A , one 
obtains: 

p + N d - n - N a = 0 (11.26) 

which is the common form of the charge neutrality equation. 

The position of the Fermi level in semiconductor : 

E c - E f = KT In for n-type(II.27)E F - E v = KT In for p-type 

VNp/ VNa/ 

(11.28) 

Ec Ec 

g 

E^E, 

Ev Ev 

intrinsic n-type 

Figure 11.25: A shift of the position of the Fermi energy in the band diagram and 
the introduction of the allowed energy level into the bandgap due to the doping. 
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Chapter III 

The SILYACO ATLAS Simulator 



111.1 Introduction 

SILVACO(SiliconValleyCorporation)is an Americancompany, headquartered inSantaClara, 
Califomia.Itis aleading provider ofbusinesssoftware strings simulationandfinite elementaided 
designforelectronics technologyTCAD(TechnologyComputer AidedDesign). These tools areused 
bycompaniesin themicroelectronicsresearch, development and designof devices. 

The edition number and the issue date is shown below 

• Edition 1 - July 1, 1993. 

• Edition 2 - March 1, 1994. 

• Edition 3 - June 1, 1994. 

• Edition 4 - October 30, 1996. 

• Edition 5 - April 30, 1997. 

• Edition 6 - November 1998. 

• Edition 7 - February 2000. 

• Edition 8- December 2002. 

• Edition 9 - September 2004. 

111.2 Silvaco Atlas Simulation Software 

Has ability toaccurately simulatea semiconductordeviceis crucialfor industry andresearch 
communities. The ATLAS simulation toolis specially designedfor 2Dand 3Dmodelingbasedon the 
physics ofsemiconductorcomponents, including electrical, optical and thermal properties.lt has 
alibrary ofpredefinedmaterialsasit allows theuser to definetheirownotherchoicematerials. 

III.2.1Atlas Inputs and Outputs 

Atlas uses two types of input file: 

•A text file that contains Atlas commands. 

•A structure file that defines the structure to be simulated. 

Atlas produces three types of output files: 
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•The runtime output that gives error and warning messages as the simulation proceeds. 
•The log file that stores voltages and currents. 

•The structure filethat stores 2D and 3D data relating to the values of solution variables 
The inputs and outputs for Silvaco Atlas are depicted in Figure III. 1 



DevEdit 



(structure and 
Mesh Editor) 



\ 



Runtime Output 



Structure Files 



> \ 



ATHENA 

(Process Simulator) 

Command File 



ATLAS 



Log Files 



/ 



(Device 

Simulator 



/ 




\ 



TonyPlot 






(Visualization 
Tool) 



DeckBuild 



(RunTime Environment) 



Figure III.l: Atlas inputs and outputs [41]. 



III.2.2Using "Silvaco- Atlas" software 

Included inthe"Atlas"physicalmodels[42]: 

• DC, AC small-signal, and full time-dependency. 

• Drift-diffusion transport models. 

• Energy balance and Hydrodynamic transport models. 

• Lattice heating and heatsinks. 

• Graded and abrupt hetero-j unctions. 

• Optoelectronic interactions with general ray tracing. 

• Amorphous and polycrystalline materials. 

• General circuit environments. 



• Stimulated emission and radiation 
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• Fermi-Dirac and Boltzmann statistics. 

• Advanced mobility models. 

• Heavy doping effects. 

• Full acceptor and donor trap dynamics 

• Ohmic, Schottky, and insulating contacts. 

• SRH, radiative, Auger, and surface recombination. 

• Impact ionization (local and non-local). 

• Floating gates. 

• Band-to-band and Fowler-Nordheim tunneling. 

• Hot carrier injection. 

• Quantum transport models 

• Thermionic emission currents. 



III.2.3Running Atlas inside Deckbuild 

To run Atlas inside Deckbuild use the following command line: 
> go atlas 



III.2.4Input file structure 

Following that command, the input file needs to follow a pattern in Figure I TT . 2 



Group 



1- Structure specification 



2. Materials model specification 



commands 

MESH 

REGION 

ELECTRODE 

DOPING 

MATERIAL 

MODELS 

CONTACT 

INTERFACE 



3. Method selection 



4. Solution specification 



5 Results analysis 



METHOD 

LOG 

SOLVE 

LOAD 

SAVE 

EXTRACT 

TONYPLOT 



Figure III.2: Atlas command groups [43]. 
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III.3Structure Specification 

The structure specification is done by defining the mesh, the region, the electrodes and the 
doping levels. 



III.3.1Mesh 

The mesh used for this thesis is two-dimensional. Therefore, only x and y parameters are 
defined. The mesh is a series of horizontal and vertical lines and spacing between them. From 
Figure III.3, the mesh statements are specified. 






\ *».« — 



-o.» 



• 8.4 



'■ 1 — 
• 208 



~ I — 
-108 



x.aesh loc=-250 spac=25 
x.nesh loc=-25 spac=2.5 
x.itesh loc=25 spac=2.5 

x. mesh loc=250 spac=25 

y. mesh loc=-2.9 spac=0.01 
y.mesh loc=-2.8 spac=0.03 
y.mesh loc--2.5 spac»0.0S 
y.mesh loc=-2 spac=0.05 
y.aesh loc=-l spac=0.1 
y.aesh loc=0 spac=0.1 



100 



— I — 
200 



Figure III.3: Atlas mesh. 



The general format to define the mesh is: 

X. MESH LOCATION =< V ALUE> SPACING=<VALUE> 

Y. MESH LOCATION=<VALUE> SPACING=<VALUE> 

For example, the x.mesh starting at -250 microns has spacing of 25 microns. That means it is 
relatively coarse. The x.mesh becomes finer between -25 and 25 microns with a spacing of 2.5 
microns. The y.mesh is similarly defined. For example, at y.mesh of -2.9 microns, the spacing is 
0.01 microns. Then at location y.mesh of -2.8 microns, the spacing changes to 0.03 microns. The 
mesh is coarser at y.mesh location of -1, when the spacing is 0.1. 

A coarse or fine mesh determines the accuracy of the simulation. A coarse mesh produces a 
faster simulation, but less accurate results. A fine mesh produces a slower simulation, but more 
accurate results. The areas that have afiner mesh, therefore, are of greatest interest in the 
simulation. 
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III.3.2 Region 

After defining the mesh, it is necessary to define the regions. The format to define the regions is 
as follows: 

REGION number=<integer><material_type> / 

<position parameters> 

From Figure III.4, the code that defines the regions is identified. There are six regions defined. 
The limits of each region are explicitly identified in the x- and y-axis. The regions must then be 
given a material. 




Figure III.4: Atlas region. 



From Figure HI. 5, the code defines the material for each region. Note that the color coding 
identifies the material. The regions have vertical and horizontal lines to mark their boundaries. 
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region mm-1 
region nnn»2 
region nun -3 
region nuu=4 
region nun=5 
region nun“6 



till: 



I r 0 a HAs 



uatcrial-InGaHAs 
note rial -Vacuun 
iiateriai-Vacu'm 
uaterlal=InGaMAs 
nacerlal=TnGaNAs 
naterial «Vacuun 



Micron* 



Figure III.5:Atlas regions with materials defined [44]. 



III.3.3Electrodes 

The next structure specification corresponds to electrodes. Typically, in this simulation the only 
electrodes defined are the anode and the cathode. However, Silvaco Atlas has a limit of 50 
electrodes that can be defined. The format to define electrodes is as follows: 

ELECTRODE NAME=<electrode namexposition_parameters> 

From Figure III. 6, the electrode statements are defined for the anode and the cathode. Note that 
the cathode is defined with gold as the material. The x and y dimensions correspond to region 6 
previously defined. Meanwhile, the anode is defined at the bottom of the cell for the entire 
xrange at y=0. 




Figure III.6: Atlas electrodes [45]. 
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III.3.4Doping 

The last aspect of structure specification that needs to be defined is doping. The format of the 
Atlas statement is as follows: 

DOPING distribution typexdopant_type> / 

<position parameters> 

From Figure III.7, the doping types and the doping levels are defined. Doping can be n-type or p- 
type. The distribution type can be uniform or Gaussian. 





0 



— r~ 



doping unicorn regional n. type conc»-2€i8 
doping uniCorn region=A n. type conc=le20 
doping unlforn cegion=5 p. type conc=de!6 
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Figure III.7: Atlas electrodes [46]. 



III.3.5Materials Model Specification 

After the structure specification, the materials model specification is next. 



III.3.5.1Material 

The format for the material statement is as follows: 

MATERIAL <localizationxmaterial_definition> 

Below are three examples of the material statement: 

MATERIAL M ATERIAL=S ilicon EG300=1.1 MUN=1200 
MATERIAL REGION=4 TAUN0=3e-7 TAUP0=2e-5 
MATERIAL NAME=base NC300=4el8 

In all examples, when MATERIAL appears first, it is considered the statement. When 
MATERIAL appears a second time in the first example, it is considered a localization parameter. 
In the second and third examples, the localization parameters are REGION and NAME, 
respectively. Various other parameters can be defined with the material statement. Examples of 
these parameters are the band gap at room temperature (EG300), electron mobility (MUN), 
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electron (TAUNO) and hole (TAUPO) recombination lifetimes, conduction band density at room 
temperature (NC300), among others. 

III.3.5.2Models 

The physical models fall into five categories: mobility, recombination, carrier statistics, 
impactionization, and tunneling. The syntax of the model statement is as follows: 

MODELS cmodel flagxgeneral parameter> / 
cmodel dependent parameters> 

The choice of model depends on the materials chosen for simulation. The example below 
activates several models. 

MODELS CONMOB FLDMOB SRH 
Table Ill.lshows somemodels. 



Table III.l: Examples of models in Silvaco 



Model 


Syntax 


Notes 


Concentration Dependent 


conmob 


It is a doping versus mobility table valid for 300K only. 


Concentration and 
Temperature Dependent 


analytic 


Caughey - Thomas formula. Tuned for 77-450K. 


Carrier-Carrier Scattering 


ccsmob 


Important when carrier concentration is high 
(e.g., forward bias power devices). 


Parallel Electric Field 
Dependence 


fldmob 


Required to model any type of velocity saturation effect. 



III.3.5.3Contact 

Contact determines the attributes of the electrode. The syntax for contact is as follows: 
CONTACT NUMBER=<n> INAME=<ename>IALL 
The following is an example of the contact statement. 

CONTACT NAME=anode current. 

III.3.5.4 Interface 

The semiconductor or insulator boundaries are determined with the interface statement. The 
syntax is as follows: 

INTERFACE [<parameters>] 

The following example shows the usage of the interface statement. 

INTERFACE X.MIN=-4 X.MAX=4Y.MIN=-0.5 Y.MAX=4 \ 

QF=lelO S.N=le4 S.P=le4 
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111.3.6 Numerical Method Selection 

After the materials model specification, the numerical method selection must be specified.There 
are various numerical methods to calculate solutions to semiconductor device problems. There 
are three types of solution techniques used in Silvaco Atlas: 

• decoupled (GUMMEL) 

• Fully coupled (NEWTON) 

• BLOCK 

The Gummel method solves for each unknown in turn, keeping the other variables constant and 
repeating the process until a stable solution is achieved. 

• This method tolerates relatively poor initial guess. 

• Generally, it is useful where the system of equations is weakly coupled but has only linear 
convergence. 

• It cannot be used with lumped elements or current boundary conditions. 

The Newton method solves the total system of unknowns together. The Block method solves 
some equations fully coupled while others are decoupled. It is useful when lattice heating or 
energy balance equations are included. 

The following example shows the use of the method statement. 

METHOD GUMMEL NEWTON 

111.3.7 Solution Specification 

After completing the numerical method selection, the solution specification is next. Solution 
specification is broken down into log, solve, load, and save statements. 

III.3.7.1 Log 

Log files (.log) store the terminal characteristics calculated by Atlas. These are current and 
voltages at electrodes in DC simulations. In transient simulations, the time is stored. In AC 
simulations, the small signal frequency, the conductances and capacitances are saved. 

The following shows an example of the LOG statement. 

LOG OUTFILE=myoutputfile.log 

The example saves the current- voltage information into myoutputfile.log. 
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III.3.7.2Solve 

The SOLVE statement follows the LOG statement. SOLVE performs a solution for one or more 
bias points. The following is an example of the SOLVE statement. 

SOLVE B 1=10 B3=5 BEAM=1 SS.PHOT SS.LIGHT=0.01 \ 

MULT.L LREQUENCY=le3 LSTEP=10 NLSTEP=6 

B1 and B3 specify the optical spot power associated with the optical beam numbers 1 and 3, 
respectively. The beam number is an integer between 1 and 10. BEAM is the beam number of 
the optical beam during AC photo-generation analysis. SS.PHOT is the small signal AC analysis. 
SS. LIGHT is the intensity of the small signal part of the optical beam during signal AC photo- 
generation analysis. MULT.F is the frequency to be multiplied by FSTEP. NFSTEPS is the 
number of times that the frequency is incremented by FSTEP. 

III.3.7.3Load and Save 

The LOAD statement enters previous solutions from files as initial guess to other bias points. 
The SAVE statement enters all node point information into an output file. 

The following are examples of LOAD and SAVE statements. 

SAVE OUTF=SOL.STR 

In this case, SOL.STR has information saved after a SOLVE statement. Then, in a different 
simulation, SOL.STR can be loaded as follows: 

LOAD INFILE=SOL.ST 

III.3.8Results Analysis 

Once a solution has been found for a semiconductor device problem, the information can be 
displayed graphically with TonyPlot. Additionally, device parameters can be extracted with the 
EXTRACT statement. To generate the Extract statement: 

Select Extract —Device. . .from the Commandsmenu 

In the example below, the EXTRACT statement obtains the current and voltage characteristics of 
a solar cell. This information is saved into the IVcurve.dat file. Then, TonyPlot plots the 
information in the IVcurve.dat file. 

EXTRACT NAME="iv" curve(v. "anode", I. "cathode") / 

OUTFILE="IVcurve.dat" 

TONYPLOTIVcurve.dat 
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Chapter IV 

Results and discussions 



IV.lIntroduction 

This work is to study the transport of electrical charges in a solar cell based on Indium Gallium 
Nitride (InGaN). By varying the alloy composition (Ga replaced by In), InGaN can reach all 
values of bandgap between 3.40 eV and 0.7eV, which covers almost the entire solar spectrum. At 
present most of the studies of InGaN photovoltaics are conducted on hetero structures. The aim 
of this study is therefore to better understand the influence of each parameter of the solar cell for 
an improved optimization of performance. In this chapter we will present the simulation results 
of the effect of the In proportion on the InGaN solar cell performance. 



IV.2 Model description 

In this work, calculations wereall performed under 1-sun AMI. 5 illumination, that is with an 
input power of Pj n = 1000 W/m'and a temperature of 300 K using one p-njunction as shown in 
figure IV. 1. 

light 

1 111 I 

^ ~~ 
p-InGaN 

n-InCJaIM 



Figure IV. 1: Structure of an InGaN solar cell. 



IV.3 Simulation parameters 



A general expression for the wavelength-dependent absorption coefficient [47] of In x Gai_ xN used 
in the simulation is given as 



2i i/2 

a(cm _1 ) = 10 5 x C(E Ph - E g ) + D(E Ph - E g ) J (IV.l) 

Where E p h is the photon energy, C and D are composition dependent parameters given in table 
IV.l for In x Gai_ x N[48]. 
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Table IV. 1: Fitting parameters used to calculate the absorption coefficient of In x Gai. x N 
alloys [49]. 



Indium composition 


C 


D 


1 


0.69642 


0.46055 


0.83 


0.66796 


0.68886 


0.69 


0.58108 


0.66902 


0.57 


0.60946 


0.62182 


0.5 


0.51672 


0.46836 


0 


3.52517 


0.65710 



Earlier, equation (1) has been used to get best fit with experimental data of absorption coefficient 
of InGaN material [49]. After thecompensationvalues C and D, E p h . We findthe following 
relationship 

a(J0 = 2.2x10 5 J 1 ' 24 /A~ E a (IV.2) 

Eg is the composition-dependent energy bandgap of In x Gai_ x N (at room temperature) which is 
calculated as follows: 



EglnGaN ^EglnN + (1 x)EgG a pj bx(l x) (IV.3) 

Where x is the indium content, Eg(InN) =0.7 eV, Eg(GaN)=3.40 eV and the bowing parameter b 
is 1.43 [50], 

Other modeling parameters such as permittivity (s/so), affinity (yjand effective masses 



(m e /mh)are calculated as follows[51]: 

e^Ga^N) = 15. 3x + 8. 9(1 - x) (IV.4) 

/(I^Ga^N) = 2. 666x 2 + 3. 659x + 4. 31 (IV.5) 

meClnxGa^xN) = 0. 12x + 0. 2(1 — x) (IV.6) 

m h (In x Ga 1 _ x N) = 0. 17x + 1. 0(1 — x) (IV.7) 



The wavelength-dependent refractive index model of Adachigiven in [51] is used for In x Gai_ x N 
as follows: 




Where A(In x Ga 1 _ x N) = 13. 55x + 9. 31(1 - x) andB(In x Ga!_ x N) = 2. 05x + 3. 03. (1 - x) 
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A general expression (as shown below) for low field mobility developed earlier by M 
Farahmandet al [51] is used for electron and hole mobility of III- nitride materials: 



Ho(N,T) — Mmin,i[T/30o] 



PI , (H max,i M-min,i ^[ T / 30 q] 



P2 



+ 



1 + 



r( T /3oo^ 4 



(IV .9) 



N ref( T / 3 00) P3 



Where i = (n,p) or electrons and holes, respectively andN is the doping concentration, N re f is 
fixed at 1 xlO 17 cm -3 , T is the temperature, (31, (32, (33, (34 and y are specific parameters for a 
given material [5]. For GaN, p max and p min values for electron (hole) mobility are 1460 (170) 
and 295 (5) cm 2 V -1 -s -1 [51, 52] respectively. Similarly, p max and p min values are 3138 (340) and 
774 (15)cm 2 V -1 .s -1 [51], respectively, for electrons (holes) for InN material. Using these data, |! 0 
(N, T) is first calculated for GaN and InN material and mobility for InGaN are then linearly 
interpolated for electrons and holes as follows: 

^(InxGai-xN) = 524x + p n (GaN) (IV.10) 

ju h (In x Ga 1 _ x N) = 6. 5x + p h (GaN) (IV.ll) 

The GaN parametres, at 300 K, used in the simulation are presentedin the table IV. 2. 



Table IV.2: The simulation parameters for GaN at 300 K 



GaN used 


Parameter Value 


Energy bandgap (ev) 


3.4 


Relative permittivity £ 


8.9 


Electron mobility (cm 2 eV-ls -1 ) 


1000 


Hole mobility (cm 2 eV -1 s -1 ) 


170 


Effective density of states in the conduction 
band Nc (cm -3 ) 


2.3xl0 18 


Effective density of states in the valence 
band Nv (cm -3 ) 


4.6xl0 19 


N A (cm -3 ) 


10 18 


N d (cm -3 ) 


10 17 


Wp(pm) 


0.01 


W„ (pm) 


3 



Table IV. 3 summarises the parameters of InGaN used in the simulation 
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TableIV.3: The simulation parameters for In x Gai_ x N materials at 300 K. 



Parameters 


Expression 


Reference 


Energy bandgap (eV) 


0.7x + 3.7(1— x) - 1.43.x.(l— x) 


[51] 


Electron mobility (cm 2 eV -1 s -1 ) 


524x + (GaN) 


[52] 


Hole mobility (cm 2 eV -1 s -1 ) 


6.5x -1- [i h (GaN) 


[52] 


Relative permittivity £ 


15.3x+8.9(l— x) 


[52] 


Electron affinity y(eV) 


3.659x+2.666x 2 +4.31 


[52] 


Effective density of states in the 


0.9x+2.3(l-x) 


[52] 


conduction band N c (cm -3 ) 






Effective density of states in the 


5.3x +1.8 (1— x) 


[52] 


valence band Nv (cm -3 ) 






Absorption coefficient a(cm -1 ) 


2.2x10^1-2%-^ 


[52] 


Carrier lifetime r(ns) 


1 ns 


[46] 


Surfacerecombination velocities 


1000 


[53] 


Sp 0 ,SpL,Sn 0 9 S„L(cmS ) 







In table IV.4,we showthe values of the InGaN band gapfor different proportions oflndium. 

Table IV.4: The InGaN band gap for different In fraction. 




IV.4 Results and discussion 

Reliable solar cell modeling should be based on accurate and reliable optical and electrical 
modeling parameters of the device. Hence a prerequisite to solar cell modeling is the use of a 
photon energy-dependent optical constant in the form of n(E) (see equation (IV. 8)) and k(E) 
values where 

k(E) = a • X/4n (IV.12) 

ais absorption coefficient given in equation (IV. 1) of the solar material. Such dispersion relations 
describing n(E) and k(E) of the III-N network for wide energy gaps are usually unavailable in the 
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commercial SILVCO TCAD software package. We have imported these values from Adachi 
[54], 

We have varied the In proportion to study its effect on the solar cell performance. A typical 
SILVACO structure output for x=0.39 is shown in Figure IV. 2. 

ATLAS 

Data tram np Jrv39&«61N Gel str 
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14 
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Ms .tee o »« aoo 

■MM 

Figure IV.2: The TCO/p- InGaN/ n-InGaN sola cell with x=0.39. 
IV.4.1 Results for proportion ofindium x=0 (GaN) 

The current density-voltage (J-V) characteristic for x=0 is shown in Figure IV. 3. 




Figure IV.3: The J-V characteristic for x=0. 
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The output power-voltagecharacteristic for the x=0 is shown in Figure IV.4. 




Figure IV.4: The P-V characteristic for x=0. 

The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 5. 



Table IV.5: Parameters of p- InGaN/ n-InGaN solar cell with x=0. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V m ax 

(V) 


P 

L max 

(mW/cm 2 ) 


FF 

(%) 


n 

(%) 


0.67 


2.76 


0.63 


2.63 


1.65 


89.3 


16.50 



forthis percentageofindium (GaN); wehavelessdensity of currentandvoltagegreater,this is due 
tothe breadth in bandgap (3.4 eV),Whichabsorbedin domain (UV). 

IV.4.2 Results for proportion ofindium x=0.20 (Ino. 2 Gao. 8 N) 

The current density- voltage (J-V) characteristic for x=0 is shown in Figure IV.5. 
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Figure IV.5: The J-V characteristic for x=0.20. 

The output power-voltagecharacteristic for the x=0.20 is shown in Figure IV. 6. 




Figure IV.6: The P-V characteristic for x=0.20. 



The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV.6. 

Table IV.6: Parameters of p- InGaN/ n-InGaN solar cell with x=0.20. 



j sc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


Vmax 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


0 

(%) 


11.19 


2.11 


10.20 


1.99 


2.03 


85.73 


20.31 
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We notea significant increasein (J sc ) and solar cell efficiency!] (The highestvalue), but 
smalldecrease in V oc . 



IV.4.3 Results for proportion ofindium x=0.39 (Ino.39Gao.6iN) 

The current density-voltage (J-V) characteristic for x=0.39 is shown in Figure IV.7. 




Figure IV.7: The J-V characteristic for x=0.39. 

The output power-voltagecharacteristic for the x=0.39 is shown in Figure IV. 8. 




Figure IV.8: The P-V characteristic for x=0.39. 



The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV.7. 
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Table IV. 7: parameters of p- InGaN/ n-InGaN solar cell with x=0.39. 



J sc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V max 

(V) 


P 

1 max 

(mW/cm 2 ) 


FF 

(%) 


n 

(%) 


12.75 


1.52 


11.59 


1.39 


1.61 


83.35 


16.11 



We observe by increasing the percentage of indium(x=0.39) a decrease in all parameters, except 
density of current(increasedless than x=0.2). 

IV.4.4 Results for proportion ofindium x=0.52 (Ino. 52 Gao. 48 N) 

The current density- voltage (J-V) characteristic for x=0.52 is shown in Figure IV. 9. 




Figure IV.9: The J-V characteristic for x=0.52. 



The output power- voltagecharacteristic for the x=0.52 is shown in Figure IV. 10. 
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The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 8. 



Table IV.8: parameters of p- InGaN/ n-InGaN solar cell with x=0.52. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V max 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


n 

(%) 


13.2 


1.15 


11.96 


1.03 


1.24 


81.38 


16.11 



Also observed that by increasing the percentage of indium; all parameters decrease exeeptJ sc . 



IV.4.5 Results for proportion ofindium x=0.57 (Ino.57Gao.43N) 

The current density-voltage (J-V) characteristic for x=0.57 is shown in Figure IV. 11. 
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Figure IV.ll: The J-V characteristic for x=57. 



The output power-voltagecharacteristic for the x=0.57 is shown in Figure IV. 12. 




Figure IV.12: The P-V characteristic for x=0.57. 



The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 9. 
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Table IV.9: parameters of p- InGaN/ n-InGaN solar cell with x=0.57. 



Jsc 

(mA/cm 2 ) 


V oc (V) 


Jmax 

(mA/cm 2 ) 


V max 

(V) 


P 

L max 

(mW/cm 2 ) 


FF 

(%) 


n 

(%) 


13.31 


1.02 


12.02 


0.92 


1.10 


80.91 


11.06 



Noteatthe beginning ofthis percentagethe stabilitycurrent (Increase isalmost non-existent (+0. 1 
mA)). V oc decreasesless thanthe previous. 

IV.4.6 Results for proportion ofindium x=0.65 (Ino.65Gao.35N) 

The current density- voltage (J-V) characteristic for x=0.65 is shown in Figure IV. 13. 




Figure IV.13: The J-V characteristic for x=0.65. 



The output power-voltagecharacteristic for the x=0.65 is shown in Figure IV. 14. 
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Figure IV.14: The P-V characteristic for x=0.65. 

The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 10. 



Table IV. 10: parameters of p- InGaN/ n-InGaN solar cell with x=0.65. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V max 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


0 

(%) 


13.43 


0.84 


12.10 


0.73 


0.89 


79.27 


8.95 



By more increasing of the percentage of indium; there is no change on the density of short- 
current(Increase isalmost non-existent (+0.1 mA)) but all other parameters have slightly 
decreased. 

IV.4.7 Results for proportion ofindium x=0.68 (Ino. 68 Gao. 32 N) 

The current density-voltage (J-V) characteristic for x=0.68 is shown in Figure IV. 15. 
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Figure IV.15: The J-V characteristic for x=0.68. 



The output power-voltagecharacteristic for the x=0.68 is shown in Figure IV. 16 




Figure IV.16: The P-V characteristic for x=0.68. 



The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 11. 

Table IV.11: parameters of p- InGaN/ n-InGaN solar cell with x=0.68. 
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Jsc 

(mA/cm 2 ) 


V oc (V) 


Jmax 

(mA/cm 2 ) 


V max 

(V) 


P 

L max 

(mW/cm 2 ) 


FF 

(%) 


n 

(%) 


13.46 


0.77 


12.06 


0.67 


0.82 


78.25 


8.20 



Increasing of the percentage of indiumestimatedbya small amount (0.03); there is also asmall 
changein the density of short-current (Increase) but all other parameters have slightly decreased. 

IV.4.8 Results for proportion ofindium x=0.78 (Ino. 78 Gao. 22 N) 

The current density- voltage (J-V) characteristic for x=0.78 is shown in Figure IV. 17. 




Figure IV.17: The J-V characteristic for x=0.78. 



The output power-voltagecharacteristic for the x=0.78 is shown in Figure IV. 18. 
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Figure IV.18: The P-V characteristic for x=0.78. 



The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 12. 



Table IV.12: parameters of p- InGaN/ n-InGaN solar cell with x=0.78. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V m ax 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


0 

(%) 


13.51 


0.57 


12.03 


0.48 


0.58 


75.72 


5.89 



Increasingof the percentage of indiumestimateby a largeamount (10%);all parameters decrease 
except J sc , that is great value (E g = 1.1 eV). 

IV.4.9 Results for proportion ofindium x=0.84 (Ino. 84 Gao.i 6 N) 

The current density- voltage (J-V) characteristic for x=0.84 is shown in Figure IV. 19. 
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Figure IV.19: The J-V characteristic for x=0.84. 



The output power-voltagecharacteristic for the x=0.84 is shown in Figure IV. 20. 




Figure IV.20: The P-V characteristic for x=0.84. 

The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 13. 



Table IV.13: parameters of p- InGaN/ n-InGaN solar cell with x=0.84. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


Jmax 

(mA/cm 2 ) 


V m ax 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


0 

(%) 


13.50 


0.46 


11.96 


0.39 


0.46 


73.61 


4.66 
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Increasingof the percentage of indium to alloy; all parameters have slightly decreased. 



IV.4.10 Results for proportion ofindium x=0.89 (Ino. 89 Gao.nN) 

The current density-voltage (J-V) characteristic for x=0.89 is shown in Figure IV. 21. 
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0,00 0,05 0,10 0,15 0,20 0,25 

Voltage (V) 



0,30 



0,35 



0,40 



Figure IV.21:The J-V characteristic for x=0.89. 

The output power-voltagecharacteristic for the x=0.89 is shown in Figure IV. 22. 




Figure IV.22: The P-V characteristic for x=0.89. 
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The extracted external parameters of this solar cell from the J-V and the P-V characteristics are 
presented in table IV. 14. 



Table IV.14: parameters of p- InGaN/ n-InGaN solar cell with x=0.89. 



Jsc 

(mA/cm 2 ) 


Voc(V) 


J max 

(mA/cm 2 ) 


V m ax 

(V) 


Pmax 

(mW/cm 2 ) 


FF 

(%) 


0 

(%) 


13.46 


0.39 


11.68 


0.32 


0.37 


71.10 


4.66 



By more increasing of the percentage of indium; all parameters have slightly decreased. 



IV.4.11 Results for proportion ofindium x=[0-0.89] 

The current density-voltage (J-V) characteristic forx= [0-0.89] is shown in Figure IV. 23. 




voltage (V) 



-■-0,89 
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-•-0 




3,0 



Figure IV.23: The J-V characteristic for different In proportion (x=0— >0.89]. 

This figure all cases presented individually before are assembled together to show the overall 
effect of the In proportion. When weincrease the percentage ofindium, there is arapid increasein 
the short circuit current density (J) anddecrease inopen circuit the voltage. 

Increasing the proportion ofindium (x)decreases the bandgapE g and thusreduces theabsorption of 
light(or in the absorption coefficient a) (see equation (IV. 1. a)), whichaffectsthe density of 
current (increasing) (J) and decrease voltage (V). 
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IV.5Summary of indium content in the InGaN (x) absorber layer effect 

IV.5 .IThe effect on the short circuit current 

The effect of theincreaseof indium (In) content in the InGaN on the short circuit current are 
summarized in Figure IV. 24. It is observed that In increases with the increases of x. This is as 
mentioned before; x/=> Eg\=> a e-h pair/=> J sc /. 
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Figure IV.24:Influenceof indium content in the InGaN layer on the density of short circuit 
current of the TCO/p- InGaN/ n-InGaN cell simulated by SILVCO. 



IV.5 .2 The effect on the open-circuit voltage V oc 

The effects increases of indium (In) content in the InGaN on the open-circuit voltage are 
summarized in Figure IV. 25. It is observed that Indecreases with in increases x. This is as 
mentioned before; x /=> Eg\=> a \=> e-h pair/=>J sc /=> V oc \. 
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Figure IV.25:Influenceof indium content in the InGaN layer on the open-circuit voltage of 
the TCO/p- InGaN/ n-InGaN cell simulated by SILVCO. 

IV.5 .3 The effect on the fill factor (FF) 

Theeffect of increases in indium (In) content in the InGaN on the fill factor is summarized in 
Figure IV. 26. It is observed that Indecreases with in increases x. 
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Figure IV.26:Influenceof indium content in the InGaN layer on the fill factor of the TCO/p- 

InGaN/ n-InGaN cell simulated by SILVCO. 
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IV.5 .4 The effect on the efficiency 

The effects of increasingthe indium (In) content in the InGaN on the efficiency are summarized 
in Figure IV. 27. It is observed that supervaluewhen x=0.2 (V 0C =2.11;J SC =11.19 mA/cm ).This 
isbecause ithas the largest power (P max ) -Compared with (GaN), which has thelargest (v oc =2.76 V), 
but JscVery small (J sc =0.67 mA/cm"), 




Figure IV.27:Influenceof indium content in the InGaN layer on the efficiency of the 
TCO/p- InGaN/ n-InGaN cell simulated by SILVCO. 

IV.6 Comparison between best results of p-n and n-p at x=0.78, In 0 . 78 Ga 0 . 22 N 

The current density-voltage (J-V) characteristic for p-n x=0.78 and n-p x=0.78 is shown in 
Figure IV. 28. 
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Figure IV.28:The J-V characteristic for a p-n and n-p junction InGaN solar cells for 

x=0.78. 



Comparison between best results of p-n and n-p at X=0.20, h^oGagoN represented in table IV. 15. 

Table IV. 15: Comparison between a p-n and n-p junction InGaN solar cells for x=0.78. 



The parameters 


p-n junction 


n-p junction 


Density of short-circuit J sc (mA/cm 2 ) 


13.51 


14.40 


Open-circuit voltage V oc (V) 


0.57 


0.61 


Fill factor FF 


75.72 


76.58 


Efficiency (%) 


5.89 


6.80 



It is observed thatJ sc ,V oc and rpFFof then-pj unction are better than of the p-n junction. 
This isdueto the mobility difference; /r n >^ /l ; (reduce the recombination carriers (e-h)). 
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General Conclusion 

In this work we have studied a p-In x Gai_ x N/n-In x Gai_ x N, homo-junction solar cell. The aim is to 
study the effect of the Indium proportion on its characteristics (the short current density J sc , the 
open circuit voltage V oc , the fill factor FF,the efficiency r| and the maximum output power 
PMAx)-The proportion ofindium for (In x Gai_ x N)was varied from 0 to 0.89. The study was carried 
out using numerical simulation under the SILVACO software.The In x Gai_ x N layer thickness is 
fixed inboth sidesat W a =0.01pm and W d =3 pm respectively. Thedopingof the two sides 
wereN a =10 18 cnT 3 and N d =10 17 cm -3 ) respectively. 

Forthe x=0 proportion the short circuit density was quite low (0.67 mA/cm") while the open 
circuit voltage has a large value (2.76 V). A fill factor was around 89.3 which give an efficiency 
of 16.5%. The addingof 20 % Indium improves the short circuit current significantly to about 
11.19 mA/cm". Thesolar cell efficiency^ also improved while asmalldecrease in V oc was 
observed (2.11 V). 

After thatwith the increase inthe percentageindiumit was noteda slight increase in J sc untilup to 
asaturated value (J sc =13.51 mA/cm') whenthe ratio of Indium is x=0.78. The other parameters all 
decrease. The In x Gai_ x N p-n junction solar cellswith different In compositions have been 
simulated. It isobserved that an In 0 . 78 Gao. 22 N(E g =l.leV) p-n junctionsolar cell exhibits the 
highest density of currentl3.51 mA/cm 2 (Voc =0.57V andq=5.89% , FF=73.61%). 

It has to be noted though that the aim is to increase in the short circuit current since this 
homojonction solar cell will be used in multijunction solar cell and in these types of solar cells 
the short circuit current is the most important parameters. 

We have also compared this structure to an n-p structure with the same parameters. It was found 
that the latter is slightly better. So as a suggestion for further work, is to study the effect of the In 
proportion of the n-p based solar cell. It is also worth to study the effect of this homojuntion 
solar cell in the multijunction solar cell. It may be also interesting to study a heterojuntion solar 
cell to find if there is any improvement. 
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